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Abstract 


We discuss a class of exact finite energy solutions to the vacuum source-free Maxwell field equations as models for multi- 
and single cycle laser pulses in classical interaction with relativistic charged test particles. These solutions are classified 
in terms of their chiral content based on their influence on particular charge configurations in space. Such solutions 
offer a computationally efficient parameterization of compact laser pulses used in laser-matter simulations and provide a 
potential means for experimentally bounding the fundamental length scale in the generalized electrodynamics of Bopp, 
Lande and Podolsky. 
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Advances in laser technology have made possible the 
exploration of physical processes on unprecedented tem¬ 
poral and spatial scales. They have also opened up new 
possibilities for accelerating charged particles using laser- 
matter interactions. Multi- and single cycle high inten¬ 
sity (10^° — 10"® Watts/cm^) laser pulses can be produced 
using Q-switching or mode-locking techniques [T]. Such 
pulses can accelerate charged particles such as electrons 
to relativistic speeds where radiation reaction and quan¬ 
tum effects may influence their dynamics. Lower intensity 
pulses have also been used as diagnostic tools for exploring 
the structure of plasmas in various states mm- In order 
to interpret experimental data involving classical laser in¬ 
teractions with both charged and neutral matter, theoret¬ 
ical models gHl] rely crucially on parameterizations of the 
electromagnetic fields in laser pulses, particularly in situa¬ 
tions where traditional formulations using monochromatic 
or paraxial-beam approximations have limitations [HHin]. 

In this Letter we discuss a viable methodology for pa¬ 
rameterizing a particular class of propagating solutions to 
the source free classical Maxwell equations in vacuo that 
offers an efficient means to explore the classical effects of 
compact laser pulses on free electrons in dynamical regimes 
where quantum effects are absent. The parameterization 
is based on a remarkable class of explicit solutions of the 
scalar wave equation found by Ziolkowski [IIHIS] follow¬ 
ing pioneering work by Brittingham m- Such solutions 
can be used to construct classical Maxwell solutions with 
bounded total electromagnetic energy and fields bounded 
in all three spatial directions. With simple analytic struc¬ 
tures their diffractive properties can be readily determined 
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together with the behaviour of charged particle-pulse inter¬ 
actions over a broad parameter range without recourse to 
expensive numerical computation. Finally, we argue that 
such parameterizations can be used to find compact finite 
energy solutions to other linear wave equations. This is 
illustrated by showing that the generalized theory of Bopp 
HZ], Lande [T5] and Podolsky m admits such particular 
solutions that reduce to the Maxwell solutions when a fun¬ 
damental length parameter in their theory tends to zero. 
Compact laser pulses in this theory might be used to ex¬ 
plore properties of the theory by searching experimentally 
for bounds on this parameter. 

If a complex scalar field a satishes □ a = 0 on space- 
time and is any covariantly constant (degree 2) anti¬ 
symmetric tensor field on spacetime (i.e. = 0) for 

all yt,iy,S = 0,1, 2, 3, then the complex tensor field = 
— di,A^ satisfies the source free Maxwell equations 
in vacuo with: 

A, = {a n^;3) JW\ (1) 

where jgl is the determinant of the spacetime metric and 
u denotes the Levi-Civita alternating symbol. In the 
following g refers to the Minkowski metric tensor field, in 
which case the components can be used to encode 
three independent Hertz vector fields and their dual^ 
General solutions to □ a = 0 can be constructed by 
Fourier analysis. In cylindrical polar Minkowski coordi¬ 
nates {t,r, z,0}, axially symmetric solutions propagating 
along the z-axis have, for z > 0, the double integral repre- 


^In the language of differential forms on Minkowski spacetime 
A = -*-^(0:11), F = dA where di^da = 0, the 2-form 11 satisfies 
Vn = 0 and ★ denotes the Hodge map associated with g. 
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Figure 1: Power profile P of the (CM, 1) laser pulse at Z = 0, T = 0 
with parameters {A = 600, 'I'l = 1, 'I '2 = 1000, $ = 0.001, H = 1} 

sentation a{t,r,z) = dw e“*“*5(a;, r, z) where: 
a{uj,r,z) = ( k /t.j(fc) dk 

Jo 

in terms of the zero order Bessel function and the speed of 
light in vacuo c. 

Conditions on the Fourier amplitudes fc^{k) can be 
given so that the Hertz procedure above gives rise to real 
singularity free electromagnetic fields with finite total elec¬ 
tromagnetic energy. A particularly simple class of pulses 
that can be generated in this way follows from the complex 
axi-symmetric scalar solution: 

£2 

+ {ipi + i{z - ct)) {'ip 2 - i{z + ct)) 

where £oi V’li i /’2 are arbitrary real positive definite param¬ 
eters with physical dimensions of length. The relative sizes 
of '01 and 02 determine both the direction of propagation 
along the z—axis of the dominant maximum of the pulse 
profile and the number of spatial cycles in its peak magni¬ 
tude. When 01 02, the dominant maximum propagates 

along the z—axis to the right. The parameter io deter¬ 
mines the magnitude of such a maximum. The structure 
of such solutions has been extensively studied in [5D1 [U] in 
conjunction with particular choices of H^;^ together with 
generalizations discussed in [251 US] • 

In general the six anti-symmetric tensors with compo¬ 
nents in a Minkowski Cartesian coordinate system 

are covariantly constant and can be used to construct a 
complex eigen-basis of antisymmetric chiral tensors 
with s G {CE, CM} and k G {—1, 0,1}, satisfying 

(3) 


where the operator represents 0 rotations about the 
z—axis generated by —ide on tensor^ These in turn 
can be used to construct a complex basis of chiral eigen- 
Maxwell tensor fields F®’”. The index s indicates that the 
CE (CM) chiral family contain electric (magnetic) fields 
that are orthogonal to the z—axis when k = 0. The chi¬ 
ral eigen-fields F®’ ° inherit the axial symmetry of a{t, r, z) 
while those with k = ±1 do not. The directions of elec¬ 
tric and magnetic fields in any of these Maxwell solutions 
depend on their location in the pulse and the concept of 
a pulse polarisation is not strictly applicable. The chiral 
content as defined here can be used in its place. Non-chiral 
pulse conhgurations can be constructed by superposition 
with arbitrary complex coefficients C®’ 
The energy, linear and angular momentum of the pulse 
in vacuo can be calculated from the components of 
the Maxwell stress-energy tensor — — 

where = Re(F^i,). If e and lo denote time- 
dependent real electric and magnetic 3-vector fields asso¬ 
ciated with any pulse solution, its total electromagnetic 
energy J', for a fixed set of parameters and any z, is cal¬ 
culated from 

J = — f dt f {ex h)-dS (4) 

PO J — oo J S 

where S can be any plane with constant z = zg > 0. For 
spatially compact pulse helds in vacuo this coincides with 
the total pulse electromagnetic energy 

p poo p‘2.'k poo 

£ = pdV= dz dO rdr p{t,r, z,6) (5) 

Jv J — oo Jo Jo 

where p = ^ ^eg e • e + is integrated over all space 
V. This follows since V • (e x b) = —po 9t p- To corre¬ 
late J with other laser pulse properties and the choice of 
parameters, we bring the pulse into classical interaction 
with one or more charged point particles. The world-line 
of a single particle, parameterized in arbitrary coordinate 
as = C^(t) with a parameter r, is taken as a solution 
of the coupled non-linear differential equations 

Mr) = ^^Mr))V^{r) ( 6 ) 

in terms of the particle charge q and rest mass mg, for 
some initial conditions ^(0),F(0), where the particle 4- 
velocity satisfies V' = —1 and its 4-acceleration is ex¬ 
pressed in terms of the Christoffel symbols F*^^^ as Afi = 
dr V^{t) + Vs{t) Vpij) ^l{i{T))■ In the following, radi¬ 
ation reaction and inter-particle forces are assumed neg¬ 
ligible. From the solution ^(r) one can determine the in¬ 
crease (or decrease) in the relativistic kinetic energy trans¬ 
ferred from the electromagnetic pulse to any particle and 
the nature of its trajectory in the laboratory frame. This 


^In terms of the Lie derivative, Oz = —H^8g s^nd = d{x± 

iy) A dt, = dzA where x = rcos{e),y = 

r sin(6) 
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Figure 2: Three-dimensional spacecurves for particles subject to an incident (CM, 1) laser pulse (left), (CM,0) laser pulse (centre) and 
(CM,—1) laser pulse (right) with parameters {A = 600, = 1000, <I> = 0.001, S = 1}. Each particle has initial velocity {R(0) = 

0, 9(0) = 0, Z(0) = shaded circular disc region indicates the initial spot size (R = 10000 for (CM, ±1) laser pulses and R = 20000 

for a (CM, 0) laser pulse) relative to the black markers on the spacecurves that denote the initial positions of the charged test particles. 


information can then be used to correlate the dynamical 
properties of the interaction with the laser pulse proper¬ 
ties fixed by the parameters. To facilitate this exercise, 
it proves important to reduce the above equations of mo¬ 
tion to dimensionless form and fix the physical dimensions 
of the fields involved. The Minkowski metric tensor field 
9 = gfiudx^ dx'' (with g^i, = diag(— 1 , 1 , 1 , 1 )) in inertial 
coordinates x° = ct,x^ = x,x‘^ = y,x^ = z) has MKS 
physical dimensions [L]^. The MKS dimension of electro¬ 
magnetic quantities follows by assigning to eoT)ji/ dx'^^dx’'^ 
in any coordinate system the physical of electric charge 
dimension. Furthermore, in terms of Minkowski polar co¬ 
ordinates {t,r,z,6}, introduce (for ease of visualization) 
the dimensionless coordinates {R = ^ ~ ~ 

and dimensionless parameters A, 'I'j = (j = 1,2) where 
[4'j] = [d>] = [S] = 1, [£o] = [L]. Then with the dimension¬ 
less complex scalar field a{T,R,Z) = a{t,r,z) and greek 
indices ranging over {T, R, Z, 9} with ^ 

write 




(7) 


for a choice of dimensionless covariantly constant tensor 
n ^/3 so that [eoAfj_ dx^] has the physical dimension of elec¬ 
tric charge and 


(.2ir 


J = 


£ = 


f 


dT dR dep{T,R,z,e) 


dZ£{T, Z). 


The parameter A controls the strength of all electric and 
magnetic fields in Fps for fixed values of the parameters 
'I' 2 , <f>, S and the overall scale io will be fixed in terms of 
the total electromagnetic energy of the pulse. For a choice 
of such parameters the real fields e and b enable one to cal¬ 
culate a numerical value F such that = ^QT. The diffrac¬ 
tion of the pulse peak along the z—axis can be used to 
define a pulse range relative to the maximum of the pulse 
peak at z = 0. To this end, the density £(T,Z) defines 
the dimensionless range Zrg by £{0,0)/£{Ti, Zrg) = 2, 


where the peak at Z = Z^g > 0 and T = Ti > 0 is 
half the height of the peak at Z = 0, T = 0. If dur¬ 
ing the interval [ 0 ,Ti] the pulse propagates with negli¬ 
gible deformation in Z, one may estimate its width Z^ 
at half height and the dimensionless pulse axial speed 
/3 = Zrg/Ti- This yields the dimensionless pulse dura¬ 
tion or temporal width Tq = Z^jfd. From these dimen¬ 
sionless values one deduces the pulse MKS characteristics 
in terms of £o and hence J. If the picosecond is used as 
a unit of time, the pulse duration becomes to = ioTojc = 
ioZ^/{f3c) = N 10“^^ sec for some value N and hence io = 
{cf3N/Z^)lO~^^ metres, J = (T(3cN/Z^) Joules, 

Zrg = io^Zrg = {'E.[3cNZrg/Z^u) 10 “^^ metres and z^ = 
'B.c/dn 10“^^ metres. A dimensionless spot-size of the pulse 
at Z = Zo > 0, T = Zo//3 is then determined by the be¬ 
haviour oiP{R, Zo/P, Zo,9). At each value of Zq this func¬ 
tion of R and 9 has a clearly defined principal maximum. If 
one associates a circle of dimensionless radius Rs{Zo) with 
such a maximum locus it can be used to define a spot-size 
with radius rs(zo) =io^Rs{Zo) = (c/3Af$i?s(Zo)/^u,)IO“^^ 
metres at z = zq. 

Figure [2 displays a clearly pronounced principle maxi¬ 
mum in the power density profile P as a function of X = 
Rcos{9) and Y — Rsin{9) at Z = 0,T = 0 for a spe¬ 
cific choice of the parameters (A, 'ki, '^> 2 , ‘h, S). The same 
parameter set is used to numerically solve ^ for a col¬ 
lection of trajectories for charged particles, each arranged 
initially around the circumference of a circle in a plane or¬ 
thogonal to the propagation axis of incident CM type laser 
pulses with different chirality. The resulting space curves 
in 3-dimensions, displayed in figure clearly exhibit the 
different responses to CM pulses with distinct chirality val¬ 
ues. The instantaneous specific relativistic kinetic energy 
of a particle with laboratory speed u is 7 — 1 in terms of 

the Lorentz factor 7 given by 7 “^ = In figure 3 

this quantity is displayed on the left for a charged partiae 
accelerated by a fixed chirality (CM,—I) type pulse where 
the pulse energy is varied by changing A. On the right the 
energy transfer dependence on pulse chirality for both CE 
and CM type pulses with fixed laser energy is displayed. 
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Figure 3: Specific kinetic energy transfer to any charged test particle. On the left, the (CM,1) pulse has parameters {^i = 1, ^^2 = 
1000, $ = 0.001, H = 1}. On the right, relative values of 7 — 1 are displayed for various (s, k) pulses with parameters {A = 1, ’Fi = 1, ^2 = 
1000,$ = 0.001, H = 1}. The differences between the energy transfers for some pulses appear indistinguishable relative to others owing to 
the logarithmic scales employed. In all cases, the charged particle has initial position {R{0) = 1,^(0) = J,Z(0) = 1} and initial velocity 
{i?(0) = 0,Z(0) = ^,0(O) = O}. 


We deduce that the momentum and angular momentum 
[24] in the propagation direction can transfer an impulsive 
force and torque respectively to charges lying in an orthog¬ 
onal plane. More generally, the classical configurations of 
a high energy pulse labelled CE and CM can be distin¬ 
guished by their interaction with different arrangements 
of charged matter. 

Furthermore, by a suitable choice of parameters, (CE, 
k) type modes can be constructed that yield the same 
physical properties , Zrg, z^, 13) for all k. Similarly the 
(CM,«;) type modes yield a k independent set with physical 
properties distinct from those determined by the (CE,k) 
modes. The pulse group speed magnitudes (as defined 
above) of all these configurations are determined numer¬ 
ically and are bounded above by the value c. To illus¬ 
trate some of these statements, table summarizes the 
MKS laser pulse characteristics for a specific choice of 
{Ti, '^ 2 , d), S} and various values of A. 

Motivated by the desire to ameliorate the divergences 
in perturbative QED a number of generalized theories of 
electromagnetism have been proposed. To date, there is 
little experimental evidence for testing their predicted de¬ 
partures from Maxwell’s theory. However, with the in¬ 
crease in laser technology one may now be entering regimes 
that may discriminate between such theories. In particu¬ 
lar, Bopp-Lande-Podolsky electrodynamics is linear in the 


s 

CE 

CE 

CM 

CM 

CM 

K 

0 

±1 

0 

±1 

±1 

S (J) 

8.07 

7.69 

8.41 

7.69 

76903 

A 

40000 

1 

0.05 

1 

100 

Zrg (m) 

1.44 

2.43 

0.899 

1.80 

1.80 

to (ps) 

3 

3 

2 

3 

3 

Zia (mm) 

0.899 

0.899 

0.600 

0.899 

0.899 

r,(0) (m) 

0.018 

0.009 

0.012 

0.009 

0.009 

I (TWcm-2) 

2.646 

1.009 

0.930 

1.009 

10090 


Table 1: Table showing MKS laser characteristics for various {s,k) 
laser pulse configurations with parameters {iPi = 1, 'I '2 = 1000, >I> = 
0.001, S = 1}. 


electromagnetic fields, contains a fundamental length Aq 
and approaches Maxwell’s theory when this length tends 
to zero. Unlike Maxwell’s theory, it contains solutions de¬ 
scribing static point charges with finite electromagnetic 
energy. 


The classical source-free Bopp-Lande-Podolsky field 
equations in vacuo for the complex electromagnetic field 
tensor = dg,A^ - are: 

= 0 . ( 8 ) 

Clearly, all classical source-free vacuum Maxwell solutions 

satisfy these equations. However, there exist additional 
solutions to ([^ that are not source-free vacuum Maxwell 
solutions. For example, such additional solutions of (|^ 
follow from Q provided: 

'Aa + X^'A'^a = 0. (9) 


Furthermore, it follows that as Aq —>■ 0, one recovers the 
class of Maxwell solutions discussed above. For non-zero 
Aq, there exist particular non-Maxwellian solutions to (|^ 
satisfying 

□ a -I- = 0. 

In Minkowski cylindrical coordinates {t,r, 0, a non- 
separable axi-symmetric compact laser solution is given 
for arbitrary real Ag by: 


a{t,r, z) 
where 



ar,t,z) \ 

■^0 J 


C{t,r, z) 


C(t,r, z) 


ipi + i{z — ct) 


■02 — i{z + ct) 


in terms of arbitrary real positive-definite parameters 0i, 
02 that determine propagation along the positive z-axis 
and where Ki{z) denotes the first order modified Bessel 
function of the second kind. The function 0(t, r, z) is posi¬ 
tive for all real r, t, z and hence both the real and imaginary 
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parts of a{t,r,z) are bounded. 

The interaction of the Bopp-Lande-Podolsky pulse with 
classical charged test particles follows from the divergence 
of the Bopp-Lande-Podolsky stress-energy-momentum ten¬ 
sor |25j . This in turn leads to the same equation of motion 
Q for the test particle trajectories [5^ and may offer an 
experimental means to discriminate between the Maxwell 
and Bopp-Lande-Podolsky descriptions of source-free elec¬ 
tromagnetic fields in vacuo. 

Further details of this approach using scalar field mod¬ 
ulated Hertz potentials based on the construction Q for 
finding spatially compact solutions to other linear tensor 
and spinor wave equations will be presented elsewhere. 
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